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ABSTRACT 

We examine the evolution of the spatial counts-in-cells distribution of galaxies and 
show that the form of the galaxy distribution function does not change significantly as 
galaxies merge and evolve. In particular, bound merging pairs follow a similar distri- 
bution to that of individual galaxies. From the adiabatic expansion of the universe we 
show how clustering, expansion and galaxy mergers affect the clustering parameter b. 
We also predict the evolution of b with respect to redshift. 

Subject headings: galaxies: statistics — cosmology: theory — large-scale structure of 
universe 



1. Introduction 



The galaxy spatial distribution function f(N, V) is a simple but powerful statistic which char- 
acterizes the locations of galaxies in space. It includes statistical information on voids and other 
underdense regions, on clusters of all shapes and sizes, on filaments, on the probability of finding 
an arbitrary number of neighbors around randomly located positions, on counts of galaxies in cells 
of arbitrary shapes and sizes randomly loc ated, and on galaxy correlation fu nctions of all orders. 
These are just some of its representations (|Saslawll2000l : ISaslaw fc Yangll2010l ). Moreover it is also 
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closel y related to the di s tribution function of t he peculiar velocities of galaxies around the Hubble 
flow (|Saslaw et al.lll99ol : Eeong fc Saslawll2004l ). 



A physically motivated form of the distribution function for galaxies in quasi-equilibriu m was 



derived and generalized using a statistical mechanical approach. (lAhmad et ah 



also been generalized to systems containing particles of two different masses (jAhmad et al 



2002 



2006b ) It ha s 



2006a|) 



In its simplest form, the probability of finding N galaxies in a cell of volume V is given by the 
gravitational quasi-equilibrium distribution(GQED): 

N(l - b) 



fv(N) 



(N(l-b) 



Nb) N 1 exp (-JV(1 - b) - Nb) 



(1) 



where N is the average number of galaxies in a cell and the clustering parameter b = —W/2K is 
the ratio of the gravitational correlation energy W to twice the kinetic energy K of peculiar veloc- 
ities relative to the Hubble flow. Althoug h other distributions have been proposed (e.g. negative 
binomial; for an early review see t hey ar e generally not physically motivated. Extensive 

computer simulations (summarized in ISaslawl 120001) designed to test the GQED ag ree closely with 



the a nalytical results. Observations at both low (jSivakoff fc Saslawll2005l ) and high (jRahmani et al 



2009 ) redshifts suggest that the form of the distribution function /y(JV) is essentially unchanged 
over a wide range of redshifts. However, mergers among galaxies could have modified the form of 
fy(N) significantly compared to the simple model which involves no mergers at all, and yet they 
did not. In this paper we examine the robustness of the GQED to merging galaxies. 



2. The positions of merging galaxy pairs 

Merging galaxy pairs are often extended structures that would have a different interaction 
potential than simple spheres or point masses. With a different interaction potential, they may 
also be distributed differently. Extended structures resulting from mergers will not only have a 
different interaction potential, they also have a different mass and should be treated differently. To 
consider how these extended structures influence the GQED, we first consider how they modify the 
interaction potential. Using the modified potential we can then derive the distribution function for 
an ensemble of species with a range of masses. 

The generalized i nteraction potentia l between particles (galaxies) each of which has an isother- 
mal halo is given by (jAhmad et al.ll2002l ) 

,/ x Gm 2 , . 

= >2 + e2) l/ 2 (2) 

where G is the gravitational constant, r is the separation between a pair of particles, m is the mass 
of each particle, and e is a parameter related to the finite size of a galaxy in proper coordinates. 
We can generalize this interaction potential by factoring out the e terms to get 

Hr) = -— «(e/r) (3) 
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where n{e/r) represents a modification to the Newtonian potential. This mo dification only affects 
the potential energy part of the configuration integral given by equation (3) of lAhmad et al.l (|2002) 



Qn{T,V) 



exp [- 



(ri,r 2 ,. 



, r„)T~ 



d 



■3AT 



(4) 



where T is the kinetic temperature of the ensemble in units where Boltzmann's constant is unity, 
a nd 4> is the interparti cle potential energy of the ensemble. By following the procedure in section 2 
of lAhmad et al.l (|2002l ). the potential energy part of the Hamiltonian for a 2-galaxy system becomes 



Q 2 (T,V) = V 2 



1 + 



3Gm 2 



2T(n)- 1 /3 



c 



e 

R~i 



(5) 



where R\ is the scale where the two-galaxy correlation function is negligible, n is the number of 
particles per unit volume, and 



c 



(± 



Ri 2r 



describes how a modification to the potential changes the partition function. 



(6) 



The mod i ficati on given by equation (jBJ) is analogous to a(e/Ri) given by equation (16) of 
Ahmad et al.l (|2002l ). but can describe a generalized modification of the potential rather than the 
particular modification that arises from an isothermal halo. The effects of a modified potential 
enter into the distribution function only through the parameter b, which now becomes 



(3/2)G 3 m%T- 3 C(e/ J Ri) 
1 + (3/2)G 3 m 6 nT- 3 C(e/ J Ri 



(7) 



For an attractive potential, n(e/r) is is always positive and hence Q{e/R\) > for all values of R\ 
so from equation (J2]), < b t < 1. 



From the thermodynamic variables of the system given by equations (26)- (30) of lAhmad et al 



(|2002l ). we see that the forms of the distribution function /y(iV) and the thermodynamic functions 
of the system are essentially unchanged by a modified potential although their values are different. 
This shows that the form of the galaxy distribution function is robust to modified potentials. This 
also allows us to treat merging galaxy pairs as a separate species and extend the GQED to describe 
the distribution of merging galaxy pairs. 

A merging galaxy pair which is bound will have its dynamical center at its center-of-mass, and 
will relax to form a single galaxy at its center-of-mass. Hence the position of the centers-of-mass 
of these merging pairs will be related to the positions of the merged galaxies. While they merge, 
these merging pairs form an extended system with an external potential given by the vector sum 
of the external potentials of both galaxies in the merging pair: 



Mr) 



GmiM Gm 2 M 



Gm 2 



r + xi 



r + x 2 



r + x i 



+ 



r + x 2 



(8) 
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where r is the distance from the center of mass of the merging pair to a more distant galaxy 
of mass M, and xi and X2 are the distances from each component to the center-of-mass of the 
system. To approximate the case when mergers between similar-sized galaxies have the dominant 
effect on the distribution function, we assume that all galaxies have the same mass m so that 
m = m\ = 77i2 = M. From equation ([8]), the first order approximation to the external potential 
contains a factor of 2m. In the case where ||r|| ^> ||x||, the first order term dominates and we 
obtain the modification factor K2 ~ 2 to the potential that arises from treating these merging pairs 
as single extended particles. 

The universe however does not contain only merging pairs. To model the presence of individua l 
galaxies that are not currently merging, we consider a two-species distribution (jAhmad et al.ll2006al ) 
where one species consists of merging pairs, and the other species consists of galaxies that are not 
merging. For simplicity we consider a system containing two species of different extended particles 
where species 1 represents individual galaxies with an average mass of m, each with a halo, and 
species 2 represents bound pairs of merging galaxies with a total mass of 2m. Using the modification 
to the GQED described above, the potential between a pair of particles of species 1 (each of which 
is a galaxy) is 

Gm\ 



4>i( r ) 



-«i(ei/r) 



0) 



where K\{e\Jr) is a softening facto r that arises from an extended halo with a physical extent 
described by e\ (e.g. isothermal halo (jAhmad et al.ll2002l )). Likewise, the external potential between 
a particle of species 1 (a galaxy) and a particle of species 2 (a bound merging pair) is given by 



Gmim 2 
h\r) = K2(e 2 /r) 



(10) 



where m,2 = 2mi because there are two galaxies of mass m\ in each merging pair, and the mod- 
ification to the potential is given by 1^2(^2/ r )- Since a modified potential only changes b in the 
single-species distribution function, a modified potential only changes the two-species clustering 
parameter b m in the two-species distribution function 



bm. — 



Ni /3nT- 3 



N 1 + (3nT- 3 + N 1 + /3i 2 nT- 3 1 + 



N 2 f3 12 nT- 3 



_b f (N 2 /N 1 )(f3 1 2/(3) \ 

N2/N1 V 1-&+OWW 



where j3 and f3\2 are given by 



P12 



RiJ 



-(Gm 
2 V 



2\3 
l) 



2r 



— } dr 



(Gm 1 m 2 fC2 



_f2_ 

Ri 



-{Gmim2) 



Ri 2r 

R\ 



2 K 2 



dr 



and b is the single-species clustering parameter given by ([Ahmad et al.ll2002l ) 

/3nT- 3 



1 + pnT 



-3 4 



(11) 

(12) 
(13) 



(14) 
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The two-species distribution function is thus (jAhmad et al.ll2006al ) 



wo - Fa -<>> 

xexp[-iV(l -b m ) -Nb m \. 



N(l -b) + ((3u/P)N b 
l-b+(/3 12 /(3)b 



i N 2 



(15) 



The two-species distribution reduces to the single-species distribution function in the limit N2/N1 — > 
and ~~ ^ 1- When Aq S> N2 and N is large, the deviation from the GQED is small because 
the [N(l — b) + Nb\ 1 term dom inates. Measurements of the fraction of merging pairs in the 



VVDS catalog (Ide Ravel et al.1 120091 ) suggest that for redshifts of z < 1, N2/N1 « 10%. Assuming 
these mergers are between galaxies that are similar in mass, 1712 ~ 2mi so P12/P ~ 8. For large 
cells which are a representative sample of the universe, A" > 100 and the difference between the 
two-species distribution and the single-species distribution is small on the level of about 5%. Hence 
under these conditions, the only significant effect of galaxy mergers in this context is a change in 
the average mass of a galaxy. 



Redshift evolution of b 



To determine the change in the clustering parameter b we consider a merging pair of galaxies 
each of mass m which approach each other with velocities of vi and V2. Since momentum is 
conserved, the merged galaxy follows the trajectory of the center of mass of the progenitors, and 
has a final velocity after the merger of Vf = (vi + V2)/2. The final velocity of the merged galaxy 
depends on the detailed dynamics of the system, but by averaging over all orientations, we find 
that mergers will not change the average kinetic energy of an ensemble. Hence the more important 
contribution to the evolution of b comes from the change in the positions of galaxies and the 
expansion of the universe. 



We extend the analysis of ISaslawl (|1992l ) to describe an ensemble where galaxies merge by 
considering the effect of the a diabatic expansion of the universe . The equations of state for internal 



energy U and pressure P are (jSaslaw 



2000 



Ahmad et al.ll2002l ) 



U = -NT(1 



2b) 



and 



NT 



(1-6) 



(16) 



(17) 



where V is the volume and T is the kinetic temperature of peculiar veloci ties in units where the 
Boltzmann constant is 1. Equations (TH|) and (fTTl) can be combined to get (ISaslawl ll992) 



tttT'— 3 



b^nT 



1 + b nT 



n T 1 3 



b PT- 



(18) 
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The analysis of section [2] generalizes equation (|18p to extended objects with the generalized form 
of b given in equation (|7|). By comparing equations (|18j) and ([7]) we see that bo is given by 

where m is the average mass of a galaxy, and £ is a function of order unity that depends weakly 
on e/R\. 

Because R\ is defined as the scale where the two-galaxy correlation function becomes negligible, 
the universe is approximately uniform averaged over scales > R±. Here we take R± to be the scale at 
which the two-galax y correlation function begins to decrease faster than a power law. Measurements 



from the 2DFGRS ( Hawkins et al.ll2003l ) have indicated that R\ is about 12/i -1 Mpc at which the 
two-galaxy correlation function is of the order 10 -2 . Here h = Ho/ WO is the reduced Hubble 
constant. We note that such cells are large enough to contain individual field galaxies and clusters 
of galaxies, and hence would be an approximately representative sample of the universe. Assuming 
that on such scales, galaxies have isotropic average velocities, then for cells with a radius larger 
than R\, galaxies are as likely to enter a cell as they are to leave a cell. With this assumption, the 
total mass in each comoving cell M c would be approximately constant, and each cell would have 
on average N = M c /m galaxies and d(rnN) = 0. Therefore N oc m _1 and bo can be written as a 
function of N instead of m. This transforms equation f|18f) into the form 



b = bo{N)PT-* = \G*(^f\ cQr)/'/ m 



Differentiating equation (|20p with respect to N gives 



b N \ 6Ri V dN J T)P J N K S ; V J 
from which we define the term 

. (22) 

S 6Ri V dN y TjP V ' 

In general, we do not rule out the possibility that (,(e/Ri) may indirectly depend on N, hence the 
<91n£(e/i?i)/<9iV factor in may be nonzero. 



In the case of adiabatic expansion, equations (fTUj) and (fT7|) . give 

3 



2 

Equation (i20"j) implies 



= dU + PdV = -(1-26) TdN + NdT\ Wp - ZNTdb + NT(1 - b) — . (23) 



dV 
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and hence using dV/V = 3dR/R where R is the scale length of the universe, we have 

= — 2b)TdN - |(1 - 2b)NT^ - 3NTdb + 3iVT(l - 6)^. (25) 
Rearranging the terms and using equation (|21|) , we find in terms of redshift z oc 1/R — 1 

* 1-b A + 66 + l-26y 1 (26) 



dz l + z \ 8b 

To illustrate how mergers contribute to the time evolution of b, we consider a simple model 
of a galaxy. In our model, galaxies have isothermal halos with a characteristic radius e, and all 
galaxies have the same density so that in a cell of total mass M c 

e = a(mf 3 = a^y 3 (27) 
for some constant of proportionality a such that e depends on th e aver age mass of a galaxy. We 



use the GQED and form of £ for such a case from I Ahmad et al.l (|2002l ) to obtain the constraint 

(28) 



-1/18 < C* < 0. The extremes of this constraint gives 

db 1-6/246 



dz 1 + z V 5 + 146 

for the case where C* = an d 

,//, 1-6/ 1 3G6 ^ (2Q) 



dz l + z V29 + 786 
for the case with = —1/18 with all other cases occurring in between equations (|28p and ([29 



Since < 6 < 1, we compare the two cases numerically and find that the difference between the 
two cases is small at less than 2%. This result tells us that although galaxy mergers can influence 
the time evolution of 6, their influence mostly results from changes of the number of galaxies in a 
cell. 



4. Conclusion and Future Work 

We have established that the effects of galaxy mergers leave the form of the galaxy distribution 
function essentially unchanged and just alter the parameters of the counts in cells distribution. In 
particular, by describing bound merging pairs as objects with a modifie d gravitational pote ntial, 



we obtain a modified form of the two-species counts in cells distribution ([Ahmad et alj|2006al ) and 
show that it only changes the counts in cells distribution slightly from the single-species result given 
by equation ([T]). 

As a result of mergers, the clustering parameter 6 increases with time, and we have shown that 
it depends very weakly on the physical extent of a galaxy and the scale R\ at which the two point 
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correlation function is negligible. The effect of the physical extent of a galaxy changes db/dz by 
less than 2%, which shows that the evolution of b depends mainly on the adiabatic expansion of 
the universe and the change in the number of galaxies from mergers. 

These results show that even when we take galaxy mergers into account, we can not only 
reproduce the GQED bu t also trace the evoluti on of the clustering parameters. However, an analysis 
of the GOODS catalog (jRahmani et al.l 120091 ) indicates a large variation between the North and 
South fields and suggests that the sample is probably too small to draw any meaningful conclusions 
about the evolution of b at high redshift. Future surveys however may provide sufficiently large 
samples at high redshifts to test our predicted evolution of b. 
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